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Generalized-ensemble simulations enable the study of complex adsorption scenarios of a coarse¬ 
grained model polymer near an attractive nanostring, representing an ultrathin nanowire. We 
perform canonical and microcanonical statistical analyses to investigate structural transitions of 
the polymer and discuss their dependence on the temperature and on model parameters such as 
effective wire thickness and attraction strength. The result is a complete hyperphase diagram of 
the polymer phases, whose locations and stability are influenced by the effective material properties 
of the nanowire and the strength of the thermal fluctuations. Major structural polymer phases in 
the adsorbed state include compact droplets attached to or wrapping around the wire, and tubelike 
conformations with triangular pattern that resemble ideal boron nanotubes. The classification of 
the transitions is performed by microcanonical inflection-point analysis. 


I. INTRODUCTION 

Although the recent advances in nanotechnology make 
it possible to study smaller systems than ever before, 
the systematic understanding of miniature structures of 
macromolecules attached to inorganic adsorbents is an 
experimental challenge P, P]. The typical approaches of 
experimental verification of theoretical predictions and 
the theoretical understanding of experimental results is 
often limited to very specific questions On the other 
hand, computer simulations are a comparatively inexpen¬ 
sive tool that allows for comprising studies, which can 
address generic problems as well. Pertaining to polymer 
adsorption at substrates, generic properties of this tran¬ 
sition include common transition pathways from the des¬ 
orbed into the different adsorption phases or vice versa 
(“wetting”, “dewetting”), but also transitions between 
structural phases under the influence of external factors 
such as temperature and solvent quality. General mate¬ 
rial properties of the adsorbent, e.g., the propensity to 
bind macromolecules and the attraction (or repulsion) 
range, affect the binding process as well. 

There have been numerous computational studies of 
this complex behavior by means of simplihed, coarse¬ 
grained models in the recent past. Yet, there are still 
more questions than answers about the nature of the ad¬ 
sorption phases that are difficult to approach. Strictly, 
most of these structural phases are not phases in the ther¬ 
modynamic sense, since the features of these phases and 
the transitions in between them are governed or at least 
substantially influenced by finite-size effects It has 
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been a fairly profound finding that the thermodynamic 
behavior of compact phases of finite polymers is rather 
irregular due to discrete crystal symmetries and surface 
effects (3-01 • 

Most of the recent generic studies aiming at unraveling 
the structure of the phase space and the adsorption dy¬ 
namics have been performed for polymers and proteins 
interacting with a per fectly flat substrate with 
or without pattern |lll - l^ . Much less literature exists 
about polymer adsorption at fluctuating surfaces such 
as membranes [l^, [25| , because of the substantially in¬ 
creased complexity of the problem. Flexible polymers 
can also optimally ada^t to curved surfaces such as cylin¬ 
drical substrates (2fil - l^ . which causes the formation of 
polymer structures different from the adsorption phases 
at flat substrates. 

Nanotubc-polymer composites are a particularly in¬ 
teresting class of hybrid systems with potential for na¬ 
notechnological applications. Besides the individual me¬ 
chanical and electronic properties of nanotubes, in par¬ 
ticular carbon nanotubes [2^ hybrid nanotube- 

polymer systems have an even broader spectrum of ad¬ 
vanced applications in technology and medicine. These 
features have been exploited, for example, in photon¬ 
ics and nanotubes coated with conducting polymers 
were used to design biosensors [s^ ■ Potential medical ap¬ 
plications include serum protein-coated gold nano-rods 
for the selective targeting of cancer cells [H, [s^. The 
understanding of the polymer wetting behavior of nano¬ 
tubes [ 35 II is a key to further developments in the techno¬ 
logical preparation of such hybrid systems. The investi¬ 
gation of the interaction of macromolecules with carbon 
nanotubes has also been the subject of numerous com¬ 
putational studies [sgI - Ij^ . 

Recently, we investigated the extreme case of poly¬ 
mers adsorbed on ultrathin nanowires and constructed 
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a structural phase diagram based on stable low-energy 
states [13 . This is particularly interesting, because 
the nanowire geometry is topologically different from 
a nanocylinder, which leads to a different monomer- 
substrate interaction model. Therefore, to consider the 
nanowire as the limiting case of a nanocylinder with zero 
radius is nontrivial. Experimentally resolved structures, 
for example for the spread of glycerol droplets on car¬ 
bon fibers and epoxy resin on aramid filaments [43 | or 
liquid droplets at thin cylinders exhibit similari¬ 

ties to the adsorbed, globular conformations we found 
in our simulations of a simple coarse-grained polymer- 
nanowire model. The outer membrane of certain bacte¬ 
ria possess an ordered, hexagonally packed surface layer 
and although the units of the outer protein layer are 
not formed by polymers, images of long tubular por¬ 
tions surrounded by hexagonally packed units [4^ re¬ 
semble polymer-tube structures we identified for strongly 
attractive nanowires [1^. Potential applications of hy¬ 
brid polymer-nanowire assemblies include tube systems 
for the transport of small molecules that can basically be 
modeled in any desired shape. 

In this paper, we extend our previous study of prop- 
erties of lowest-energy polymer adsorption phases (43| 
by introducing the temperature as the external con¬ 
trol parameter for thermal fluctuations. By means of 
generalized-ensemble Monte Carlo methods and by sys¬ 
tematic variation of temperature and model parameters 
such as effective wire thickness (volume exclusion) and 
attraction strength, we investigate the thermodynamic 
properties of the structural phases of adsorbed and des¬ 
orbed polymers and the transitions between them. The 
competition between energetic ordering and entropic di¬ 
versity that causes structural transitions will be investi¬ 
gated by the conventional canonical statistical analysis of 
thermodynamic response quantities and more advanced 
approaches such as the microcanonical inflection-point 
analysis ili3- The latter method has proven to be par¬ 
ticularly useful, whenever finite-size effects dominate or 
substantially influence structural transitions. 

The paper is organized as follows. After the introduc¬ 
tion of the hybrid polymer-wire model in Sect. im we 
briefly review the conformational phase diagram of poly¬ 
mer ground-state structures adsorbed to a nanowire in 
Sect, nni Results for the adsorption thermodynamics of 
the polymer are analyzed and discussed in Sect. IIVI The 
microcanonical analysis of the adsorption transition is 
performed in Sect.|Vl Our main conclusions are summa¬ 
rized in Sect. ED 


II. THE POLYMER-WIRE MODEL 

We employ a generic, coarse-grained model for a self¬ 
interacting homopolymer with N monomers in the prox¬ 
imity of an attractive one-dimensional stringlike sub¬ 
strate. The energy of a conformation X of the poly¬ 
mer consists of contributions for the interaction between 


nonbonded monomers, bending stiffness of the polymer 
chain, and the interaction of each monomer with the 
nanowire: 

N-2 N N-1 
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The pairwise interaction between nonbonded monomers 
is described by a standard Lennard-Jones potential. 
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where is the distance between monomers i and j. 
The potential minimum Vlj (?'““) = — Cm is located at 
r™™ = 2^/®crni. In our simulations, we used energy and 
length units for which Cm = 1 and tTm = 1, respectively. 
Bonds between monomers adjacent along the chain are 
stiff and of unit length, ra+i = 1. Covalent bond vectors 
connected to the ith monomer form a bending angle 0^. 
The bending energy 

Wend(cos6»i) = k(1 - cos6»i) (3) 


is a remnant of the AB model, which has initially been in¬ 
troduced as a coarse-grained hydrophobic-polar peptide 
model [i^ that has proven to be useful for the quali¬ 
tative description of generic features of conformational 
transitions associated with protein folding [i^ and ag¬ 
gregation processes . In this model, the bending stiff¬ 
ness is set to K = 1/4. The homopolymer version of the 
model has also been used in recent adsorption studies of 
polymers [l^ HI, [^, li^, IsH . 

For the determination of the string potential, we con¬ 
sider a monomer interacting with each infinitesimal ele¬ 
ment of the string, which is supposed to extend infinitely 
into ±z direction. We assume that this interaction is of 
van der Waals type and can be described by a Lennard- 
Jones potential. By using the cylindrical coordinate rep¬ 
resentation, the total potential of a monomer i in the 
field of the string is then given by [i^ [5l| 
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where we have introduced the perpendicular distance rx;i 
of the monomer from the string axis and the van der 
Waals radius or effective “thickness” of the string erf, 
which is related to the potential minimum distance r™™ 
via: 
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FIG. 1. Cylinder potential [Eq. ([6|] for different radii R € 
(0 ... 1] and the string potential [Eq. Q] for Uf = ef = 1 
for constant total “charge”, i.e., rjcyi = (main figure) 

and for constant “charge” density r?cvi = 1 (inset). The inset 
figure corresponds to Fig. 4 in Ref. 

The string “charge” density rjf compensates for the addi¬ 
tional dimension by the extra CTf factor that remains after 
integration over the standard form of the Lennard-Jones 
potential. It can be conveniently used to set the energy 
scale of the potential. We adjust it in such a way that the 
minimum monomer-string energy is t4tring(''’™'”) = — Cf, 
independently of CTf. In this case, rjf « 0.528/CTf. 

It is instructive to view the string potential (| 31 ) from a 
different perspective, namely as the R ^ 0 limit of the 
corresponding potential of an attractive cylinder. This is 
nontrivial, because of the different topological dimensions 
of string and cylinder. The cylinder-monomer potential 
can be written as (^ : 

p27t pR 

V'cyi(i?,r_L;i) = 7ref77cyi / d4> I pdp 

JO JO 
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where the integration in z direction has already been per¬ 
formed. The distance between the ith monomer and any 
point in the cylinder disk in-plane with the monomer is 
given by r' = {r^.^ + p^ — 2pr±-i p and cj) are 

the polar coordinates in the plane of cylinder disk and 
monomer. Compared to Ref. |26l |. we have introduced in 
the potential ([6]) the volume “charge” density rjcyi of the 
cylinder. If this density is considered a constant and is 
thus independent of R (as in Ref. [1^), the string po¬ 
tential (|4|) cannot be obtained as the i? —>■ 0 limit of the 
cylinder potential. Only if the total “charge” is fixed, 
i.e., rjcyi scales with the inverse area of the cylinder disk 
I/ttR?, then 

I)5tring(^_L;2) — lim Wyl(R, rj_.f ) . (7) 

R —^0 

For both cases, the i? —)■ 0 limits of the cylinder poten¬ 
tial are plotted in Fig. [TJ The systematic discussion of 


polymer adsorption at cylindrical substrates will be done 
elsewhere [^ . 

In our model of the polymer-wire system, the polymer 
is not grafted to the string and can move freely. The 
string, pointing into z direction, is located in the center 
of the periodic simulation box with edge lengths 2N in x 
and y directions to prevent the polymer from escaping. 

In our simulations, we performed different runs u sing 
multicanonical [52 h^ and Wang-Landau sampling [^, 
[5^ and found the results to be coherent. Conformational 
updates included local crankshaft and slithering-snake 
moves, as well as global spherical-cap and transla¬ 
tion moves M- The crankshaft move corresponds to the 
rotation of a single monomer about the axis given by 
the vector between its nearest neighbors and changes the 
conformation locally. For the slithering-snake update, a 
monomer is cut at one end of the chain and pasted at the 
other end, inverting the bond vector. This move proves 
useful to avoid trapping in very dense conformations. 
The spherical-cap update consists of a pivot rotation of a 
single bond and the subsequent shift of the monomers at 
either end of the polymer, keeping all bond lengths fixed. 
This update allows for larger steps in the conformational 
space. Finally, the global translation update allows for a 
displacement of the entire polymer chain relative to the 
string. Update types were chosen randomly with equal 
weight, which allowed for the efficient sampling of the 
conformational space. 

The chain lengths studied varied from = 30 to 
N = 100 for the analysis of the thermodynamic adsorp¬ 
tion behavior and up to = 200 for the investigation of 
ground-state properties. 


III. LOWEST-ENERGY CONEORMATIONS OE 
ADSORBED POLYMERS 

The investigation of conformational properties of 
lowest-energy states is generally beneficial for the iden¬ 
tification of compact crystalline or amorphous structural 
phases. The major ground-state morphologies of the 
polymer interacting with an attractive nanowire were 
found to be string-attached globular droplets and tube¬ 
like barrels wrapped around the string j43l |. In the first 
case, depending on the attraction strength and effective 
thickness of the wire, the string penetrates the globule 
(Gi: globular, included) or it is only loosely attached to 
it (Ge: globular, excluded). 

Most interesting is the formation of tubelike, or bar¬ 
rellike (B), polymer structures, which resemble nanotube 
morphologies [57|. This occurs if the string field is more 
attractive for most monomers than the compact crys¬ 
talline assembly and the effectively excluded string di¬ 
ameter is small enough. 

Topological transitions, common to all adsorption 
problems and typically accompanied by layering effects, 
can also be observed in the barrel phase. At very low 
CTf values, where monomers are attached so close to the 
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CTf = 0.5 I cTf = 1.5 CTf = 2.0 



FIG. 2. Typical low-energy conformations of a polymer with 
N = 200 monomers adsorbed at a string with effective thick¬ 
nesses CTf-values 0.5 and 1, respectively (first column), 1.5 
(second) and 2 (third) for attraction strengths et = 1, 2,3, 4, 
and 5 (top to bottom). Monomers are marked with different 
shades according to their distance from the string. 

string that they can feel the attractive field of monomers 
of the opposite side of the string as well, the polymer 
“tube” looks rather like a one-dimensional chain with zig¬ 
zag (or sawtooth) pattern. In the other extreme case, 
if the effective string thickness is too large, additional 
monomers will attach to the two-dimensional monolayer 
tube, essentially forming a new layer. Three-dimensional 
assemblies with two or more layers are also known from 
atomic nanotubes [s^. 

For comparatively large string thickness and suffi¬ 
ciently high adsorption strength, the polymer wraps only 
partially around the string. This structural phase is 
very typical for polymer adsorption transitions at cylin¬ 
ders [23 [13 and has been called “clamshell” phase (C). 

Figure [5] shows typical conformations in these dif¬ 
ferent structural low-energy phases for a polymer with 


200 monomers at various values of string parameters erf 
and Cf, supporting our general findings for other chain 
lengths |43l |. Since the monomer-string interaction is 
attractive, the polymer is adsorbed to the string in all 
low-energy phases. 

In the following, we will investigate the influence of 
thermal fluctuations upon the the formation of compact 
adsorbed polymer structures and the thermodynamic 
properties of the adsorption transition. 


IV. THERMODYNAMICS OF POLYMER 
ADSORPTION AT STRINGLIKE SUBSTRATES 

For the understanding of the thermodynamic behav¬ 
ior of the system, it is useful to perform canonical and 
microcanonical analyses of various thermodynamic quan¬ 
tities with the goal of identifying generic features of 
polymer adsorption at a stringlike nanowire for vari¬ 
ous parametrizations in the space of effective van der 
Waals radius Uf of the wire and its attractive adsorption 
strength Cf. 

Figure [3] shows heat-capacity curves Cv{T) for a poly¬ 
mer consisting of Y = 30 monomers interacting with var¬ 
ious nanowires, distinguished by their different material 
parameters CTf and ei. Most pronounced in all cases is the 
significant peak at high temperatures, which corresponds 
to the adsorption transition. For fixed values of Uf, the 
peak position, which we denote by in the following, 
scales almost perfectly linearly with the wire attraction 
strength, TX oc £{. This has already been found in 
adsorption studies at planar substrates [H, [13 and can 
intuitively be understood by the larger thermal energy 
~ fceT’ needed to release the polymer off the substrate at 
larger adsorption strengths. However, it is worth men¬ 
tioning as well that the proportionality factor depends 
nonlinearly on the effective string thickness CTf. 

This can also be seen in Fig. |4l where we have plotted 
the adsorption temperature as a function of CTf and ei for 
a 100-mer. We have also included the zero-temperature 
structural phase diagram from Ref. (43 into the T = 0 
plane to get an impression of how the compact adsorption 
phases at given CTf and Cf values look like. Not surpris¬ 
ingly, the adsorption transition temperatures depend on 
N and are different for the 100-mer, compared to the 
30-mer. However, qualitatively the generic results are 
comparable. 

After passing the adsorption transition, for tempera¬ 
tures < T < dominant polymer conforma¬ 

tions are either disordered and expanded or locally or¬ 
dered and globular. For the adsorption of the 30-mer 
at the nanowire, these phases are difficult to distinguish. 
More striking is the crystallization transition at 
into the compact adsorption phases that have been de¬ 
scribed qualitatively in Section uni As it is obvious from 
Fig. [31 the low-temperature peak positions in the spe¬ 
cific heat curves near indicating this transition, 

do not vary much and, therefore, hardly depend on the 






5 



T 


FIG. 3. Temperature dependence of the heat capacities of 
a 30-mer for different string attraction strengths (cf = 1, 3, 
and 5, respectively) and decreasing van der Waals radii: (a) 
CTf = 2.0, (b) 1.5, (c) 1.0, and (d) 0.5. Note the different 
ordinate scales. Error bars were estimated using the jackknife 
method 0]. 

material parameters of the nanowire, ef and df. Effec¬ 
tively, during this transition, the number of contacts be¬ 
tween monomers and wire does not change much and the 
way in which the polymer binds to the wire has already 
been established in the “liquid” (expanded or globular) 
phase. Note that, for strong attraction (large Cf values), 
the crystallization transition becomes weaker, the smaller 
the effective wire thickness is. This confirms the previ¬ 
ous statement that the basic, tubelike structures have al¬ 
ready formed in the “liquid” phase, with most monomers 
already in contact to the wire. 

If the wire attraction and thickness are comparatively 
small (cf < 2 and df < 1), the chain behaves basically like 
a free, flexible polymer, with minimal recognition of the 
wire, which is enclosed inside the polymer conformation 
for energetic reasons. In this case, the adsorbed liquid 
phase is preempted by the compact, crystalline phase. 
Adsorption and crystallization merge. Noteworthy are 
the very weak peaks in Fig. [3] at about Tg,™ « 1.5, which 



FIG. 4. Adsorption temperatures for different values of ef and 
effective wire thickness df for a polymer of length N = 100. 
The structural phase diagram for ground-states from Ref. 
is projected upon the T — 0 plane. 


indicate the 0 transition of the polymer chain. The weak¬ 
ness of the 0 transition for very short flexible and semi- 
flexible chains has already been observed previously 

Figure [5] shows an example for a complete structural 
phase diagram at df = 1.0 for N = 30. Representative 
conformations in the different phases at parameter values 
indicated by the circles in the phase diagram are also de¬ 
picted. The phase diagram prominently features the al¬ 
most linear adsorption transition line between desorbed 
and adsorbed phases. Transitions from the desorbed ex¬ 
tended (DE) and desorbed globular (DG) phases, which 
are separated by the well-known 0 transition line, into 
adsorbed phases lead to a variety of structures. However, 
given the rather small size of the polymer, these confor¬ 
mations can hardly be distinguished. Although we expect 
that the dominant DE structures will also adsorb as ex¬ 
tended conformations and, similarly, DG structures cross 
over into adsorbed globules, a quantitative description is 
more difficult than in the adsorption scenarios at planar 
substrates [l^ [l^ . A significant signal for the separation 
of extended and globular in the adsorbed regime could 
not be identified. The “crystallization” transition into 
compact phases dominated by the ground-state struc¬ 
tures, however, is clearly visible. In complete analogy 
to the previous analysis of ground-state structures for 
CTf = 1 in the case of the 100-mer (see also the 
CTf = 1 cross section in the projection in Fig. a , we 
clearly identify compact globular droplets with the string 
excluded (Ge) at rather small values of the string attrac¬ 
tion strength Cf. With Cf increasing, the string becomes 
included into the globule (Gi). Only for sufficiently large 
string attraction (cf > 3), tubelike “barrel” (B) struc¬ 
tures form. Further increasing Cf, these become dom¬ 
inated by monolayer (or single-walled) polymer tubes. 
The detailed discussion of the ground-state dominated 
compact phases can be found in Ref. [43 |. 
























6 






(1.0; 2.5) 


( 1 . 0 ; 1 . 0 ) 


(2.0; 1.5) 


(3.0; 2.5) 


(3.0; 1.0) 


(4.0; 1.5) '^(4.0;0.75) (5.0; 0.75) (5.0; 2.0) 


FIG. 5. Example of the structural phase diagram and repre¬ 
sentative conformations in the different conformational phases 
for = 30 at wire thickness at = 1.0. Transition tem¬ 
peratures Ttrans Were identified from peaks in the specific 
heat. The low-energy “clamshell” phase C does not exist for 

(Tf = 1.0. 


V. THE ADSORPTION TRANSITION FROM 
THE MICROCANONICAL PERSPECTIVE 


The canonical statistical analysis of different response 
quantities, including the specific heat, as functions of the 
canonical temperature is typically ambiguous. Although 
signals such as extremal points and inflection points are 
helpful indicators of thermal activity, they do not allow 
for a unique identification of transition points [ 3 , 113 , . 

Figure [5] shows the adsorption peaks in the heat-capacity 
curves for various £f values and wire thickness tTf = 1.0 
for a 30-mer and a 100-mer. 

An alternative method that has proven to be extremely 
useful is based on the fact that temperature can also be 
defined as a system property. This method is called mi- 
crocanonical inflection-point analysis and rests upon 
the shoulders of microcanonical statistical analysis . 



FIG. 6. Adsorption transition signals in the heat-capacity 
curves of (a) a 30-mer and (b) a 100-mer for different values 
of the wire attraction strength Cf and fixed effective thickness 
CTf = 1.0. For visualization purposes, only the adsorption 
peaks are shown for cf < 5.0. 


The central quantity is the density of states g{E') = 
J VX S{E' — i5(X)), which sets entropy and energy into 
relation with each other: S{E) = k-Q In g{E). It is easy to 
show that microcanonical equilibrium between two sys¬ 
tems is achieved if dSi{Ei)/dEi = dS 2 {E 2 )/dE 2 . With 
this, the inverse microcanonical temperature can be de¬ 
fined as fi{E) = dS{E)/dE. Figures (iKa) and (b) show 
for the 30-mer the microcanonical entropy and inverse 
temperature curves, S{E) and /3{E), respectively, for 
wire material parameters £f = 1.0,1.5,2.0,...,5.0 and 
(Tf = 1.0. For comparison, the same plots are shown in 
Fig. [5] for a 100-mer. 

The I3{E) curves in Figs. (iKb) and [HKb) exhibit a 
characteristic behavior that help locate transition points 
uniquely. Specific monotonic properties can even be 
used to introduce a systematic classification scheme. In 
microcanonical inflection-point analysis [47l |. the tran¬ 
sition points are defined by means of the derivative 
7 (A) = d(3{E)/dE in the following way. If 7 (A) ex¬ 
hibits a local maximum and the extremal value is posi¬ 
tive, 7 max(Atr) > 0 , then the associated inflection point 
of the microcanonical temperature at Etr must lie in a 
low-sensitivity region, in which an energetic increase has 
a minimal effect on the change of temperature within a 
certain energetic range. For a finite system, this entails 
the “backbending” of the temperature curve, which levels 
out in the thermodynamic limit and renders the Maxwell 
construction exact. The energy that spans the transition 
range is associated with a nonzero latent heat AQ > 0, 
which is why this transition is referred to as a first-order 
transition. 
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e = E/N 

FIG. 7. Microcanonical analysis of structural transitions of 
the polymer with 77 = 30 monomers interacting with wires 
of effective thickness (Tf = 1 and attraction strengths tf = 
1.0,1.5, 2.0,..., 5.0. (a) Microcanonical entropy S(e), (b) in¬ 
verse microcanonical temperature /3(e) = {l/N){dS{e)/de), 
and (c) inflection-point indicator 7 (e) = {1/N){dl3{e)/de). 


Consequently, a negative local maximum 7max(-E’tr) < 
0 , which also signals a low-sensitivity region of the /?(£’) 
curve, marks a second-order transition, because it corre¬ 
sponds to a transition point in energy space and, there¬ 
fore, AQ = 0. Note that according to these definitions 
the order of transition depends on the system size. It is 
not unusual that a transition with hrst-order signature 
in a finite system turns into a continuous, second-order 
phase transition in the thermodynamic limit. Not sur¬ 
prisingly, the order of a transition is also sensitive to 
changes of system (or model) parameters. 

Most strikingly, as can be seen in Figs.[7Kb) and (c) for 
the 30-mer, the adsorption transition that we classify as a 
second-order transition at £f = 1.0, turns into hrst order 



e = E/N 

FIG. 8 . Same as Fig. [71 but for N = 100 and attraction 
strengths ei = 1.0, 2.0,..., 5.0. 


for larger attraction strengths. The first-order charac¬ 
teristic becomes stronger with increasing £f values. The 
horizontal lines in Fig. [ZKb) represent the (inverse) tran¬ 
sition temperatures identified from the corresponding in¬ 
flection points. They coincide very nicely with the peak 
values of the specific heat curves shown in Fig. [51 The 
scenario is similar for the 100-mer. The inverse micro- 
canonical temperatures, plotted in Fig. |5Db) for various 
£f values, exhibit features of a second-order transition for 
£f < 3.0, before turning to first order at higher £f values. 
This is very similar to the previously discussed case of 
polymer adsorption at flat substrates 

The reason for the increasing entropic suppression of 
states in the phase separation regime is that less states 
are available for each energetic state in the linearly in¬ 
creasing energy interval (the energy scale grows like 
E oc £f). Since the desorption phase is almost entirely 
governed by translational entropy (in a constant energy 









































interval above the adsorption transition), and the ener¬ 
getically controlled adsorption states are shifted to lower 
energies (below the adsorption transition) with in¬ 
creasing, the energetic space in between is entropically 
“emptied” (note that the total number of states is in¬ 
dependent of Ef). For the same reason, the slope of 
the double-tangents (Gibbs construction) at the respec¬ 
tive transition points decreases for larger values of Ef, 
which effectively shifts the adsorption transition toward 
higher temperatures, which is consistent with the canoni¬ 
cal analysis of Fig.[6l Since the energy scale of adsorption 
is linear, the adsorption transition temperature scales al¬ 
most linearly as well over a large range of the adsorp¬ 
tion strength parameter, Tads ^ £f, which is a common 
feature of all minimally (additively) coupled adsorption 
models [isl lisl. 1^. 

The increased Ef range of second-order behavior of the 
100-mer, compared to the 30-mer, is an empirical indi¬ 
cation that the adsorption transition of a polymer at a 
linelike substrate is a second-order phase transition in the 
thermodynamic limit. This is consistent with the gen¬ 
eral assumption that adsorption transitions of infinitely 
large systems are continuous [^. This means that the 
width of the phase-separation regime, which is nonzero 
in the finite system, shrinks with increasing number of 
monomers and the per-particle latent heat diminishes for 
TV oo [ 13 . 

VI. SUMMARY 

We have investigated the thermodynamic behavior of 
structural phases for a hybrid model of a polymer inter¬ 
acting with stringlike substrate that resembles a nano¬ 
wire. This study extends the previous study of ground- 
state properties of the same system [13 by analyzing 
the influence of thermal fluctuations upon the adsorp¬ 
tion transition and the formation of compact adsorption 
phases. Results were obtained by means of advanced 
generalized-ensemble Monte Carlo simulations. We iden¬ 
tified all major structural thermodynamic phases from in¬ 
dications of enhanced thermal activity in the specific heat 
curves and constructed the structural hyperphase dia¬ 


gram, parametrized by temperature and string attraction 
strength. As expected, for increasing attraction strength, 
the low-temperature structures change from compact 
crystalline shapes that attach to the string or include 
it to barrellike conformations that resemble single-walled 
nanotubes. The transition into less ordered adsorption 
phases turns out to be independent of the adsorption 
strength, i.e., the transition temperature is virtually con¬ 
stant. In contrast, the adsorption/desorption transition 
temperature scales linearly with the substrate attrac¬ 
tion strength. A microcanonical inflection-point analy¬ 
sis revealed that the adsorption transition possesses an 
increasingly strong first-order characteristic, the larger 
the attraction strength is. However, it is important to 
note that this first-order behavior is a finite-system fea¬ 
ture that is supposed to disappear in the thermodynamic 
limit. It is only due to an effective, relative suppression 
of microstates in the phase-separation transition region, 
which is a finite-size effect. This does not contradict the 
common view that the adsorption transition is a second- 
order phase transition in the thermodynamic limit. 

The results presented in this paper provide general 
background knowledge needed for systematic approaches 
to the technological design of functional materials used 
for the transport of molecules (based on polymer tubes) 
and also give insight into binding options of biomacro¬ 
molecules at nanofibers (such as actin filaments) in cells. 

ACKNOWLEDGMENTS 

This work has been supported partially by the NSF 
under Grant No. DMR-I207437, by CNPq (National 
Gouncil for Scientific and Technological Development, 
Brazil) under Grant No. 402091/2012-4, and by the 
DFG (German Research Foundation) under Grant No. 
SFB/TRR 102 (Project B04). The authors also acknowl¬ 
edge support by the Jiilich/Aachen/Haifa Umbrella pro¬ 
gram under Grants No. SIM6 and No. HPC_2. Gomputer 
time was provided by the Forschungszentrum Jiilich un¬ 
der Project No. jiff39 and No. jiff43, and by the Georgia 
Advanced Computing Resource Center (GACRC) at the 
University of Georgia. Assigned: LA-UR-14-27917. 


[1] S. R. Whaley, D. S. English, E. L. Hu, P. F. Barbara, A. 
M. Belcher, Nature 405, 665 (2000). 

[2] J. J. Gray, Curr. Opin. Struct. Biol. 14, 110 (2004). 

[3] M. Bachmann, K. Goede, A. G. Beck-Sickinger, M. 
Grundmann, A. Irback, and W. Janke, Angew. Ghem. 
Int. Ed. 49, 9530 (2010) 

[4] M. Bachmann, Thermodynamics and Statistical Mechan¬ 
ics of Macromolecular Systems, (Gambridge University 
Press, Gambridge, 2014). 

[5] T. Vogel, M. Bachmann, and W. Janke, Phys. Rev. E 
76, 061803 (2007). 

[6] S. Schnabel, T. Vogel, M. Bachmann, and W. Janke, 


Chem. Phys. Lett. 476, 201 (2009). 

[7] S. Schnabel, M. Bachmann, and W. Janke, J. Ghem. 
Phys. 131, 124904 (2009). 

[8] T. Bogner, A. Degenhard, and F. Schmid, Phys. Rev. 
Lett. 93, 268108 (2004). 

[9] A. Swetnam and M. P. Allen, Phys. Rev. E 85, 062901 

( 2012 ). 

[10] M. Moddel, W. Janke, and M. Bachmann, Phys. Rev. 
Lett. 112, 148303 (2014). 

[11] T. Vrbova and S. G. Whittington, J. Phys. A 29, 6253 
(1996); J. Phys. A 31, 3989 (1998); T. Vrbova and K. 
Prochazka, J. Phys. A 32, 5469 (1999). 



9 


[12] Y. Singh, D. Giri, and S. Kumar, J. Pliys. A 34, L67 
(2001); R. Rajesh, D. Dhar, D. Giri, S. Kumar, and Y. 
Singh, Phys. Rev. E 65, 056124 (2002). 

[13] M. S. Gauso, J. Chem. Phys. 117, 6789 (2002). 

[14] J. Krawczyk, T. Prellberg, A. L. Owczarek, and A. Rech- 
nitzer, Europhys. Lett. 70, 726 (2005). 

[15] M. Bachmann and W. Janke, Phys. Rev. Lett. 
95, 058102 (2005); Phys. Rev. E 73, 041802 (2006); 
Phys. Rev. E 73, 020901(R) (2006). 

[16] N. Kallrot and P. Linse, Macromolecules 40, 4669 (2007). 

[17] J. Luettmer-Strathmann, F. Rampf, W. Paul, and K. 
Binder, J. Ghem. Phys. 128, 064903 (2008). 

[18] M. Moddel, M. Bachmann, and W. Janke, J. Phys. 
Chem. B 113, 3314 (2009). 

[19] L. Wang, T. Chen, X. Lin, Y. Liu, H. Liang, J. Chem. 
Phys. 131, 244902 (2009). 

[20] A. D. Swetnam and M. P. Allen, Phys. Chem. Chem. 
Phys. 11, 2046 (2009). 

[21] V. A. Ivanov, J. A. Martemyanova, M. Muller, W. Paul, 
and K. Binder, J. Phys. Chem. B 113, 3653 (2009). 

[22] M. Moddel, W. Janke, and M. Bachmann, Phys. Chem. 
Chem. Phys. 12, 11548 (2010); Macromolecules 44, 9013 
( 2011 ). 

[23] Y. W. Li, T. Wiist, and D. P. Landau, Phys. Rev. E 87, 
012706 (2013). 

[24] H. Popova and A. Milchev, J. Chem. Phys. 127, 194903 
(2007); ibid. 129, 215103 (2008). 

[25] S. Karalus, W. Janke, and M. Bachmann, Phys. Rev. E 
84, 031803 (2011). 

[26] A. Milchev and K. Binder, J. Chem. Phys. 117, 6852 

( 2002 ). 

[27] 1. Gurevitch and S. Srebnik, Ghem. Phys. Lett. 444, 96 
(2007); J. Chem. Phys. 128, 144901 (2008). 

[28] J. Gross, T. Vogel, and M. Bachmann, Structural Phases 
of Adsorption for Flexible Polymers on Nanocylinder 
Surfaces, to be published (2015). 

[29] M. S. Dresselhaus, G. Dresselhaus, and P. Avouris, 
eds., Carbon Nanotubes: Synthesis, Structure, Proper¬ 
ties, and Applications, Topics in Applied Physics, Vol. 80 
(Springer, Heidelberg, 2001). 

[30] A. Jorio, G. Dresselhaus, and M. S. Dresselhaus, eds., 
Carbon Nanotubes: Advanced Topics in the Synthesis, 
Structure, Properties, and Applieations, Topics in Ap¬ 
plied Physics, Vol. Ill (Springer, Berlin, Heidelberg, 
2008). 

[31] T. Hasan, Z. Sun, F. Wang, F. Bonaccorso, P. H. Tan, 
A. G. Rozhin, and A. C. Ferrari, Adv. Mater. 21, 3874 
(2009). 

[32] M. Gao, L. Dai, and G. G. Wallace, Electroanalysis 15, 
1089 (2003). 

[33] L. Wang, Y. Liu, W. Li, X. Jiang, Y. Ji, X. Wu, L. Xu, 
Y. Qiu, K. Zhao, T. Wei, Y. Li, Y. Zhao, and C. Ghen, 
Nano Lett. 11 , 772 (2011). 


[34] S. H. Brewer, W. R. Glomm, M. G. Johnson, M. K. Knag, 
and S. Franzen, Langmuir 21, 9303 (2005). 

[35] M. Q. Tran, J. T. Cabral, M. S. P. Shaffer, and A. Bis¬ 
marck, Nano Lett. 8 , 2744 (2008). 

[36] C. Wei, Nano Lett. 6 , 1627 (2006). 

[37] S. S. Tallury and M. A. Pasquinelli, J. Phys. Chem. B 
114, 4122 (2010); ibid. 9349 (2010). 

[38] W. Liu, C.-L. Yang, Y.-T. Zhu, and M.-S. Wang, J. Phys. 
Chem. C 112, 1803 (2008). 

[39] C. Caddeo, C. Melis, L. Colombo, and A. Mattoni, 
J. Phys. Chem. C 114, 21109 (2010). 

[40] E. J. Wallace and M. S. P. Sansom, Nano Lett. 7, 1923 
(2007). 

[41] E. J. Wallace, R. S. G. D’Rozario, B. M. Sanchez, and 
M. S. P. Sansom, Nanoscale 2, 967 (2010). 

[42] P. Angelikopoulos and H. Bock, Langmuir 26, 899 (2010). 

[43] T. Vogel and M. Bachmann, Phys. Rev. Lett. 
104, 198302 (2010); Phys. Proc. 4, 161 (2010) 

[44] H. D. Wagner, J. Appl. Phys. 67, 1352 (1990). 

[45] B. J. Garroll, Langmuir 2, 248 (1986). 

[46] R. Glaeser, W. Chiu, and D. Grano, J. Ultrastruct. Res. 
66 , 235 (1979). 

[47] S. Schnabel, D. T. Seaton, D. P. Landau, and M. Bach¬ 
mann, Phys. Rev. E, 84, 011127 (2011). 

[48] F. H. Stillinger, T. Head-Gordon, and G. L. Hirshfeld, 
Phys. Rev. E 48, 1469 (1993); F. H. Stillinger and T. 
Head-Gordon, Phys. Rev. E 52, 2872 (1995). 

[49] M. Bachmann, H. Arkin, and W. Janke, Phys. Rev. E 
71, 031906 (2005). 

[50] G. Junghans, M. Bachmann, and W. Janke, Phys. Rev. 
Lett. 97, 218103 (2006); J. Ghem. Phys. 128, 085103 
(2008). 

[51] T. Vogel and M. Bachmann, Comp. Phys. Comm. 
182, 1928 (2011). 

[52] B. A. Berg and T. Neuhaus, Phys. Lett. B 267, 249 
(1991); Phys. Rev. Lett. 68 , 9 (1992). 

[53] W. Janke, Physica A 254, 164 (1998). 

[54] M. Bachmann, Phys. Scr. 87, 058504 (2013). 

[55] F. Wang and D. P. Landau, Phys. Rev. Lett. 86 , 2050 
(2001); Phys. Rev. E 64, 056101 (2001). 

[56] C. Zhou and R. N. Bhatt, Phys. Rev. E 72, 025701(R) 
(2005). 

[57] T. Vogel, T. Mutat, J. Adler, and M. Bachmann, Com- 
mun. Comp. Phys. 13, 1245 (2013); Phys. Proc. 15, 87 
( 2011 ). 

[58] For a review, see C. Shen, A. H. Brozena, Y. Wang, 
Nanoscale 3, 503 (2011). 

[59] M. Bachmann and W. Janke, J. Chem. Phys. 120, 6779 
(2004). 

[60] D. H. E. Gross, Microcanonical Thermodynamics (World 
Scientific, Singapore, 2001). 

[61] E. Eisenriegler, K. Kremer, and K. Binder, J. Chem. 
Phys. 77, 6296 (1982). 


